Constraints on the cosmological relativistic energy density 
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We discuss bounds on the cosmological relativistic energy density as a function of redshift, re- 
viewing the big bang nucleosynthesis and cosmic microwave background bounds, updating bounds 
from large scale structure, and introducing a new bound from the magnitude-redshift relation for 
Type la supernovae. We conclude that the standard and well-motivated assumption that relativistic 
energy is negligible during recent epochs is not necessitated by extant data. We then demonstrate 
the utility of these bounds by constraining the mass and lifetime of a hypothetical massive big bang 
relic particle. 
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I. INTRODUCTION 



Cosmological constraints on the relativistic energy 
density (RED) of the Universe, particularly the bound 
derived from the successful prediction of the primor- 
dial abundances of D, *He, and ^Li, provide perhaps 
the strongest connections between cosmology and par- 
ticle physics. In this paper we examine constraints on 
the RED (To be specific, we consider constraints on the 
energy density of any particle that obeys an equation of 
state P = p/3.) at various cosmological epochs and dis- 
cuss how these bounds can be used to limit the mass and 
lifetime of an unstable big bang relic particle. We re- 
view the classic big bang nucleosynthesis (BBN) bound 
and the bound from the cosmic microwave background 
(CMB) anisotropy power spectrum. We present a new 
bound derived from the magnitude-redshift relation of 
Type la supernovae (SNIa), and we update the bound 
from large scale structure (LSS). Although we will at- 
tempt to derive bounds that are independent of addi- 
tional cosmological parameters, where necessary we will 
borrow constraints from other techniques in order to 
further limit the relativistic energy density of the Uni- 
verse. For example, when needed we adopt the bound 
0.1 < r^M < 0.5 (r^M is the matter density relative to 
the critical density today, pcrit = STJq/SttG), consistent 
with recent measurements of the contemporary matter 
density from lensing Q|, the rich cluster baryon fraction 

, mass-to-light ratio estimates ||^ , and the power spec- 
trum of the Lya forest 0. We also adopt a hard limit 
on the Hubble parameter {h = _ffo/100 kms^^Mpc^^) of 
0.56 < h < 0.88 in accordance with the 2a range of the 
HST Key Project 0. 
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II. CONSTRAINTS ON RELATIVISTIC 
ENERGY 



A. Nucleosynthesis Constraint 

It has long been known that an increase in the RED of 
the Universe during BBN leads to an overproduction of 
''He § (BBN reviews include 0). Prior to freeze-out of 
the weak interactions, the neutron-to-proton ratio {n/p) 
tracks its equilibrium value. As the Universe expands and 
cools, the weak interactions freeze-out and {n/p) freezes 
in. Freeze-out occurs when the weak interaction rates 
become comparable to the expansion rate H, given by 
H'^ = 8ttGp/3, where p is the mean cosmological energy 
density. Once deuterium can be synthesized nearly all 
available neutrons are incorporated into ^He. Increas- 
ing the energy density of the Universe during BBN in- 
creases the expansion rate, causing the weak interactions 
to freeze-out at a higher temperature and therefore, at a 
higher value of (n/p). An increased expansion rate also 
allows less time for free neutron decay. Both effects make 
more neutrons available for ^He synthesis. 

In addition, BBN production of ''He slowly increases 
with the baryon density or equivalently, the baryon-to- 
photon ratio, rj = us/n^. Since BBN production of deu- 
terium is a rapidly decreasing function of rj, the baryon 
density can be deduced from the D/H ratio observed in 
QSO absorption fine systems. Furthermore, because ^He 
has been produced in stars since BBN, the observed ''He 
abundance places an upper limit on its primordial abun- 
dance. Armed with this upper bound to primordial '*He 
and an independent measure of "q, we can place a bound 
on the RED during the BBN epoch {z ^ 10^). 

We illustrate the current status of the BBN bound 
to the RED as follows (see the recent analyses of |^ 
for details). We adopt a liberal range for the primor- 
dial D/H ratio that accommodates recent observations 
of deuterium in several QSO absorption line systems 
1.75 < 10^(D/H) < 5.0. We then bound the RED by 
taking a conservative bound on the "^He mass fraction, 
Yp < 0.25. This bound is consistent with both the low 
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FIG. 1: The BBN bound on the RED. The region bounded 
by the dashed curves delimits the parameter space for which 
the theoretically predicted and observed values of D/H are 
consistent. Likewise, the data on Yp are consistent with pre- 
dictions in the region below the solid curve. The intersection, 
marked by the circle, shows the maximum value of N^^^f / that 
can be made consistent with these limits on Yp and D/H. 



|pO| and high |Tl| zero-metallicity extrapolations of low- 
metallicity HII region data (but see [^). We follow 
convention by expressing limits on relativistic energy in 
terms of the energy carried by an effective number of light 
neutrino species, N^^^ff- The RED at some redshift z, 
after e+e" annihilation, is related to N^^^ff by 



eiz) 



[l + 0.135(7V,,e//-3)], 



(1) 



where we scale the RED by Pji,{z) which is the standard 
model RED carried by photons and three massless neu- 
trino species (with a current CMB temperature of 2.725 
K^,n^, = p^,{z = 0)/pcrit = 4.18/1-2 X 10-5). In 
Figure |l], we show regions of the ij-N^^eff plane allowed 
by the observed ^He and D abundances. The data re- 
quire that N^,eff^ 3.4 or, in terms of the RED during 
nucleosynthesis 6'bbn = 0{z = 10^), 



e 



BBNr 



1.05. 



(2) 



In Figure ^, we show the amount of extra relativistic 
energy permitted during BBN. 

There are two important caveats regarding the BBN 
RED bound. First, although (||) is a rather stringent 
bound on the RED, corresponding to an increase of 5% 
over the standard p^y^ it only applies during the epoch of 
nucleosynthesis. Relativistic energy injected after BBN 
is not subject to this bound (in Section III, we discuss a 
toy model in which relativistic energy is injected into the 
Universe after the light nuclide abundances have been 
fixed). Second, the BBN RED bound is flavor depen- 
dent in the sense that extra relativistic energy in the 



FIG. 2: Constraints on the RED. The horizontal lines denote 
the maximum value of Q{z) allowed by BBN, CMB, LSS and 
SNIa. Notice that the BBN constraint in the lower right-hand 
corner is hardly noticeable (i.e., very strict). By extending 
the limits to higher redshift, we have made the assumption 
that after BBN relativistic energy is only injected and never 
removed. The curved line shows Q{z) in a scenario where 
a massive big bang relic with Mf^^y^ = 1-2 x 10"^ decays 
with lifetime Tyr = 2 x lO''. There are additional CMB and 
LSS constraints on relativistic energy injected after z ~ 1000; 
however, they do not rule out the decay shown above (see 
sections III and ^ ) . 



form of degenerate electron neutrinos changes the rates 
of the weak interactions that inter-convert protons and 
neutrons thereby changing theoretical BBN yields. In- 
troducing an appropriate amount of electron neutrino 
degeneracy can always compensate for excess relativistic 
energy present during nucleosynthesis p^ , p^ , nullifying 
the BBN RED constraint. 



B. CMB Constraint 

The CMB anisotropy power spectrum is also sensitive 
to the RED of the Universe at recombination, primarily 
through the early Integrated Sachs- Wolfe (ISW) effect. 
The erosion of gravitational potentials due to incomplete 
matter domination at recombination leads to a boost in 
power, particularly near the first acoustic peak (for a re- 
view see jl^). CMB constraints on relativistic energy 
complement those from BBN in two ways. First, the 
CMB constraints are flavor independent; the power spec- 
trum measures relativistic energy regardless of its form. 
Second, the CMB constrains the RED at later epochs 
so that the CMB can be used to study the injection of 
relativistic energy after the epoch of nucleosynthesis, for 
example, by massive neutrino decays. 

Wang, Tegmark & Zaldarriaga ||l^ have compiled a 
combined CMB data set including the recent results from 
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the BOOMERANG (Ts), MAXIMA (T| and DASI ||] 
experiments. Hannestad has performed a Ukehhood 
analysis on this combined data set and found, with weak 
priors on the Hubble and tilt parameters and the baryon 
density, N^^^ff < 19 with 95% confidence. In terms of 
the RED present at recombination, 6'cmb = 0{z « 1100), 
this bound is 



the above limit to N^^eff < 17.5 at 95% confidence 
(see also [|2| ) . Kneller et al. , |^ have explored the prior 
dependence of these bounds in detail, showing that the 
bound on N^^^ff scales with the upper bound chosen for 
r^M (i.e., interesting CMB bounds on the RED require a 
prior constraint on JIm). 



0CMB< 3.2 (3) 

as depicted in Figure 2. 

A cautionary note is in order. Although the priors 
chosen by Hannestad are quite conservative, it must be 
borne in mind that the bound in does depend upon 
these priors. For example, by adopting tighter, yet rea- 
sonable, prior constraints on h and rj, Hannestad shrinks 



C. Constraints from Type la Supernovae 

Given a distant population of so-called "standard can- 
dles" , the magnitude-redshift relation is a powerful way 
to determine cosmological parameters directly |2^. In a 
standard FRW cosmology, the apparent bolometric mag- 
nitude, m(z), of a standard candle is related to its abso- 
lute bolometric magnitude, M, and redshift by 



m{z) = M + 5 log(Pi(z, l^M, f^A, f^R,)) - 5 log(i7o) + 25, (4) 

where Hq is the Hubble parameter in kms^^Mpc^^. 2?^ is the "Hubble-constant-free" luminosity distance in Mpc, 
given by 



dz 



V(l + ^)^(1 + ^Mz) + (1 + z)2(2 + z)z1]r - z(2 + z)f7A 



where c is in kms ^, f2k = 1 — ^^m — f^A — f^R, and 

E(a;) = sin(x) if ilk < 
— X if f^k = 
= sinh(x) if Jlk > 0. 

The matter, radiation and cosmological constant energy 
densities enter T>l with different powers of z, making it 
possible to utilize observations of standard candles over a 
range of redshifts to determine the cosmological density 
parameters (for a review, see 

Two groups, the Supernova Cosmology Project ]2^ 
and the High-z Supernova Search Team ||26| , have been 
engaged in a systematic study of the magnitude-redshift 
relation of high-redshift, type la supernovae in an effort 
to constrain SIm and f2A. We have re-analyzed the data 
of Perlmutter et al. |2^ allowing for a non-negligible 
contribution from JIr. For simplicity, we have assumed 
flatness to be a robust result of CMB measurements 
p?! and have performed a maximum-likelihood analy- 
sis (we take C oc subject to the priors J^m > 0, 
f^R ^ 0) in order to derive constraints on J^m and 
with = 1 — — f^R^. 



In Figure |3| we show the constraints on Q,m and Q,^ 
from SNIa. The projection f2p> = is consistent with 
earlier analyses that found 17m ~ 0.3 and J^a ~ 0.7 as- 
suming the RED to be negligible. Allowing for relativistic 
energy density, we find a degenerate set of SIm and J^r 
that are consistent with the SNIa data: the high-matter- 
content (i.e., 30% matter, 70% cosmological constant) 
flat Universe and the high-radiation-content (i.e., 20% 
radiation, 80% cosmological constant) flat Universe are 
equally good fits. Observe that regardless of whether rel- 
ativistic energy is allowed or not, the SNIa data require 
a large cosmological constant (0.6^ l^A^ 0.9). 

The 95% confidence contour is approximately fit by the 
region 

r2R + 0.6217M w 0.17 ±0.08. (5) 

It is not surprising that the data pick out this degenerate 
valley in the Qm-^^r plane because most of the SNIa data 
are from z ~ 0.5 and our degenerate valley represents 
the parameters with approximately constant luminosity 
distance at this redshift. The degeneracy can be further 
understood with the help of Figures ^ & ^. Dashed and 



^ Relaxing flatness yields a less restrictive bound on but much 
of the favored region of parameter space would be ruled out by 



age considerations, CMB measurements, or the aforementioned 
bounds on Hm- 
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FIG. 3: 68% (filled) and 95% (open) confidence contours in 
the r^M-^^R plane. For illustration, we show several labeled 
isochrones (dotted) and the elliptical 95% confidence contour 
obtained by adding an additional prior constraint of Q,m ~ 
0.3 ±0.1 (solid, heavy). Also shown are contours of constant 
effective ms for SNIa at redshifts of z = 0.5 (dashed at ms — 
23, 23.2) and z = 1.5 (dash-dot at ms = 26, 26.2). 
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FIG. 4: The high-redshift portion of the Hubble diagram for 
the SCP data (with la error bars) alongside the magnitude- 
redshift relation in several flat cosmologies: f^M = 0.3 and 
fia = (heavy, solid), flu = and S7r = 0.2 (dashed), 
Qa = 1 (dotted), JIm ~ 1 (dash-dot), Or = 1 (light, solid). 



dot-dashed lines in Figure ^ depict contours of constant 
apparent magnitude at redshifts z = 0.5 and z = 1.5 
respectively. As most of the SNIa data lie near z ~ 0.5, 
the confidence region is nearly parallel to lines of constant 
ms at z = 0.5. At higher redshifts, lines of constant 
apparent magnitude have a more shallow slope and are 
closer together, thus observations of SNIa at z > 0.5 can 
break the matter-radiation degeneracy. 

Figure ^ further illustrates the large lever arm of high- 
z SNIa for cosmological parameter estimation. Notice 
that the magnitude-redshift relation in a high-matter- 
content Universe (f^M — 0.3, H.^ = 0) is quite similar 
to that in a high-radiation-content Universe (rjpi = 0.2, 
= 0) at z^ 0.5. At higher redshift, one begins 
to probe the epoch prior to matter- A and/or radiation- 
A equality; the cosmological constant becomes increas- 
ingly unimportant compared to radiation and/or mat- 
ter and the two curves begin to diverge. The SCP 
data only extend to z = 0.83, therefore they cannot be 
used to scrutinize this earlier phase and they cannot dis- 
tinguish between a high-matter-content Universe and a 
high-radiation-content Universe. The proposed Super- 
nova Acceleration Probe ( pittp: / / snap.lbl.g"ov | ) may have 
the ability to break this degeneracy by observing many 
more SNIa at significantly higher redshift. 

Marginalizing over J^m by integrating the likelihood, 
we obtain a bound on the RED at low redshift ^sNia = 



0.5), 



OsNU^ iAh^ X 10^ (68%) or 4.8/i^ x 10^^ (95%). (6) 

In terms of the RED today, the bound is JIr < 0.14 
(68%) and f^R < 0.20 (95%). 

The SNIa constraint on the RED applies during re- 
cent epochs, namely, z^ 0.5, as can be seen in Figure 
^. Of course, in light of other estimates of cosmological 
parameters, the allowed region in Fig. ^ is not equally 
probable. Estimates of the contemporary matter den- 
sity of the Universe favor the region Q.m — 0.3 ± 0.1 
and would slightly reduce the SNIa upper bound on the 
RED. This is illustrated in Figure ^ where we also plot 
the 95% confidence contour obtained with the additional 
prior constraint U,m = 0.3 ± 0.1. This prior constraint 
on the matter density reduces the 95% bound on Q,n by 
about 25%. 

Note also that the SNIa bound on the contemporary 
RED is more stringent than bounds that follow from the 
requirement that the expansion age of the Universe be 
at least as large as the ages of the oldest objects in the 
Universe. In Figure |3[ the entire region delineated by 
SNIa data corresponds to an acceptable age. In partic- 
ular, a high-radiation-content Universe with f^M — 0.04, 
r^R = 0.20, r^A ~ 0.76 and a Hubble parameter at the 
extreme lower limit, h — 0.56, is 13.8 Gyr old. Even with 
h — 0.72, to — 10-7 Gyr, a value that is not grossly in- 
consistent with the ages of the oldest stars and globular 
clusters lEsll. 
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D. Large Scale Structure Constraints 

The effect of additional relativistic energy on the 
growth of large scale structure (LSS) has been studied 
by numerous authors |9[ |3l|, g ||. Typically, 
the introduction of "hot dark matter" was considered a 
means to suppress power on small scales and thus rec- 
oncile an = 1, Cold Dark Matter (CDM) cosmology 
with the observed power spectrum derived from galaxy 
surveys |34| . Conversely, too much relativistic energy ad- 
versely affects the growth of structure and therefore LSS 
can be used to constrain the RED. In light of mounting 
evidence, the CDM paradigm has given way to the so- 
called ACDM paradigm with ^ 0.3 and JIa ~ 0.7. 
With this in mind, we revise previous work in order to 
constrain the relativistic energy content of the Universe 
using LSS. 

The most striking feature of a CDM or ACDM 
type power spectrum is a break in the power law at, 
roughly speaking, the co-moving horizon scale at matter- 
radiation equality 



Aeq ^ Winuh)-^ h-^ Mpc. 



(7) 



This feature arises because the growth of sub-horizon 
sized perturbations is quelled by the cosmological expan- 
sion during radiation domination. Hence, perturbations 
on scales smaller than Aeq are suppressed by a factor 
~ (A/Aeq)^ relative to scales that were super-horizon 
sized at matter radiation equality. 

If the RED of the Universe is contained entirely in pho- 
tons and three light neutrino species, and if the primor- 
dial power spectrum is nearly scale-invariant, the ACDM 
matter power spectrum can be expressed in terms of only 
one quantity, the shape parameter F ~ ilMh |l35|| (this 
neglects the effect of baryons, see [Q). Several authors 
have used LSS observations on linear scales to infer ac- 
ceptable values of the shape parameter . One of the 
more permissive of these determinations is the 95% range 



Taking the lower bound F > 0.06 we immediately come 
upon a generic constraint on the mean RED: 



< 



0.06 



(11) 



With our conservative assumptions that < 0.5 and 
h < 0.88, the corresponding restriction on the RED dur- 
ing (and prior to) the epoch of matter-radiation equality, 
dhss = 0(z = 2eq), is 



dLSS^ 54 



(12) 



where 1 -I- zeq = ^m/^r- The relative weakness of 
this bound is due to our conservative lower bound on 
F. Taking the 95% band of Eisenstein & Zaldarriagga 
0' 0-15 < r < 0-58, results in 9. 

Note that this constraint ( [l^ ) allows the epoch of 
matter-radiation equality to be at redshifts as low as 
^EQ ~ 120. We can obtain a more stringent bound on 
the RED by following an argument invoked by Turner, 
Steigman & Krauss Steigman & Turner |30 and 

Turner & White Assuming a nearly scale- invariant 
primordial power spectrum, data from the COBE DMR 
experiment pof indicate that the rms density contrast 
at horizon crossing is on the order of Su ~ few x 10~^ 
pof . Meanwhile, measurements of the galaxy correla- 
tion function reveal nonlinear clustering on scales smaller 
than a critical scale, Anl bh^^ Mpc We adopt 

the conservative constraint that perturbations on scales 
smaller than Anl must have grown by at least a fac- 
tor of 7inin = 10"^ in order for the rms perturbation on 
these scales to be nonlinear. In linear perturbation the- 
ory, density fluctuations grow as S oc (1-1-2:)"^ during 
matter domination and only logarithmically during ra- 
diation domination. This implies that the matter dom- 
inated epoch must span at least three orders of magni- 
tude in redshift. With Qm < 0.5 this, in turn, imposes 
the limitation 



0.06 < F < 0.46, 



(8) 



quoted by Efstathiou & Moody , which we will use to 
constrain the RED. 

In the presence of excess relativistic energy, the horizon 
scale at matter-radiation equality is no longer given by 
(0). Rather, 



e^ssh^^ 12 or 6 



LSS 



< 



38, 



(13) 



where we have taken h > 0.56. This constraint is shown 
in Figure ^ where we summarize the constraints on rela- 
tivistic energy imposed by BBN, the CMB, LSS and the 
SNIa magnitude-redshift relation. 



Aeq ^ 16{nMh)-^^nR/n^Ji~^ Mpc (9) 
and the effective shape parameter is therefore given by 

F ~ nuhe-^^^. (10) 



III. CONSTRAINING RELIC DECAYS 

In the preceding section we reported limits on the RED 
at various epochs. In the absence of electron neutrino de- 
generacy, the BBN constraint on the RED is, by far, the 
most stringent; it leaves little room for any non-standard 
relativistic energy during the BBN epoch. One way to 
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circumvent the BBN constraint is to inject relativistic en- 
ergy after nucleosynthesis has ended, for instance, by the 
decay of a particle that was non-relativistic during BBN 
(M ^ 1 MeV) into relativistic products. As such, the 
study of relic particle decays comes part and parcel with 
constraints on relativistic energy. We examine the simple 
case of a massive, unstable big bang relic which may be 
pertinent to physics "beyond the Standard Model" and 
discuss constraints on the relic's mass and lifetime that 
follow from the bounds on the RED in Section II. It is 
important to note that, aside from the CMB constraint, 
the assumption that the relic is very massive is not crit- 
ical. In general, the particle must be non-relativistic at 
decay in order for its decay product's energy density to 
be comparable to or greater than p^^. 

Consider the decay of a massive relic particle X, with 
lifetime t, into relativistic products. Had X not decayed, 
the energy density in these particles today would be 



(14) 



In (|r^, MkeV is the mass of the particle in keV and 
Y is the ratio of the number density of the particle to 
the entropy density {Y ~ 0.039 for a light neutrino, 
y ~ 2 X 10"^" for a 5 GeV, Dirac neutrino). Given a spe- 
cific particle physics model in which X is produced in the 
early Universe, Y is fixed (in the absence of subsequent 
entropy production); however, in order to make the con- 
straints on heavy particle decays as generic as possible, 
we have chosen to keep Y as an explicit, free parame- 
ter. Further, we assume that the daughter particles are 
weakly-interacting (constraints on radiative decays are 
quite severe ^^). In all analytic calculations we assume 
that decays occur simultaneously at t = t. Energy con- 
servation during decay demands that the present energy 
density in decay products be J^d = f^x/(l + -^d), where 
ZD is the redshift at decay (i.e., z at t = r). Assuming 
the Universe to be X-dominated prior to decay. 



5.1 X 10-^Mtit^Y^^^T^/^ (15) 



where Tyr is the lifetime of X in years. 



A. CMB Constraints on Relic Properties 




FIG. 5: The regions of parameter space for a massive relic 
decay that are excluded by CMB (heavily shaded), LSS (mod- 
erately shaded) and SNIa (lightly shaded) arguments. Each 
excluded region extends to arbitrarily high MkevY. 



MlyY\yr< 7.5 



X 10" 



(16) 



For pertinent lifetimes, 10"'^^ Tyr"^ 10^, the excluded re- 
gion of parameter space is displayed in Figure ||. 

Post-recombination X decays modify the CMB power 
spectrum through the ISW effect and through shifts in 
the multipole positions of the acoustic peaks due to a 
change in the angular diameter distance to the surface 
of last-scattering. It is possible to use this modification 
of the observed CMB anisotropy power spectrum to con- 
strain post-recombination X decays but these con- 
straints would not be generic because the gravitational 
dynamics of the relic as well as its decay scheme can con- 
tribute to the ISW effect. Such constraints would have 
to be developed on a case-by-case basis considering the 
relic mass and relic abundance separately; however, we 
mention a specific case that evades the CMB bound and 
contributes a significant RED at the present epoch in see- 
In the following subsection, we show that the 



IV 



tion 

growth of LSS can place severe, yet generic, constraints 
on post-recombination decays. 



One constraint on relic properties follows from the re- 
quirement that the total energy density in relativistic 
particles JIr = fi-y^ -|- Qd, fall under the CMB bound 
in ^ during the epoch of recombination. If the particle 
decays prior to recombination the CMB constraint (|^), 
along with the analytic approximation of (plsj), implies 
that^ 



B. LSS Constraints on Relic Decays 

Large scale structure considerations also lead to 
bounds on decaying relics. Utilizing (|l5|), 9 can be ex- 



^ The CMB constraint on relativistic energy is strict enough that 
the assumption of X-domination prior to decay is untenable, 



making it necessary to integrate the equations governing heavy 
WIMP decay in a flat, Friedmann cosmology. We have performed 
the necessary integration and find that the above analytic bound 
is typically accurate to within 20%. 
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pressed as 61 ~ 1 + 1 2 . 2 Af f F'^/'^ r^/^ . Consequently, the 
effective shape parameter (O) may be written 

r ^ 17m/.[1 + \2.2MtJ:tvy^l\ll\^l^ (17) 



from which we conclude that 



(18) 



Again, taking VIm < 0.5, < 0.88 and T > 0.06, we 
deduce that 



7 ■ 



(1 + zx) 



(1 + ^eq) 



2.4 X lO^f^M/i^. 



(22) 



With our aforementioned limits on Qm and h, these scales 
grow by a maximum of 7 ~ 9.3 x 10'^. Perturbations that 
enter the horizon after X domination begins, grow by a 
smaller factor: 



7(A > Ax) ^ 9.3 X lO^f ^ 



(23) 



M, 



keV 



Y^Tyr< 16. 



Observe that (|19|) restricts the 
M^^yY'^Tyr- As has been underscored 
& Peacock and Bharadwaj & Sethi |Q 
spectrum will also exhibit a feature on 



(19) 

combination 
by McNally 
, the power 
small scales 



corresponding to the transition between radiation domi- 
nation and an early matter dominated (domination by X 
particles) phase. As the co-moving horizon scale during 
the first epoch of matter-radiation equality is given by 



Ax 



1 



17.1MkevY 



Mpc, 



(20) 



this small-scale feature can, in principle, constrain the 
combination MkevY alone. In practice, however, much of 
the interesting region of relic parameter space would cor- 
respond to nonlinear scales and such a constraint would 
require a better theoretical handle on nonlinear cluster- 
ing and bias. 

We can strengthen the bound in ( p^ by requiring that 
structure grow sufRciently. As Steigman and Turner ]3C| ] 
have noted, there are two ways in which this can occur. 
One way was mentioned in Section [I E , namely, that the 



relic decay early enough so that the most recent epoch of 
matter domination began at a redshift (1 -I- zeq) > 10"^. 
The redshift of equality is (1-|-zeq) — Om/^^d and using 
( |r^ ) and we find that this scenario requires 



Mi,vY^ryr< 0.66 



(21) 



and is relevant for lifetimes in the range 10^''^ Tyr^ 10^. 

In the presence of a massive unstable relic however, 
there can be two phases of matter domination, an early X 
dominated phase and a second matter dominated phase 
after relic decay. It may be possible for perturbations on 
scales A < Anl to take advantage of both of these periods 
of growth and thereby grow by a factor greater than 7inin- 
X domination begins at redshift (1 4- zx) = flx/^fiy, 
followed by decay at redshift (1 + zq) = flx/^D- Thus 
the total growth factor for scales smaller than the horizon 
scale at X domination. Ax, is 



Compelling 7(Anl) to be greater than 7inin forces 



MkevY< 3.6h X 10"^< 3.2 x 10"^. (24) 



Combining (|T|), 1^) and (|24|), we summarize the LSS 
constraints on early relic decays as 



M^^yY'^Tyr< 0.66 or 

Mi^yY^Tyr< 16 and MkevY< 3.2 x IQ-^. (25) 

The excluded region is shown in Figure ^. 

The above constraints (|25| ) are pertinent for lifetimes in 
the range 10"^^ Tyr^ 10^ because we assumed that the 
turnover in the power spectrum is indicative of Aeq, the 
horizon scale at the epoch of matter-radiation equality af- 
ter X decay. Alternatively, the X particles may have very 
large lifetimes, Tyr S> 10^, in which case the turnover in 
the power spectrum would be indicative of Ax, the hori- 
zon scale at the first epoch of matter-radiation equality, 
prior to X decay. In this case, the effective shape param- 
eter is simply 



F ~ (17x + f^M)/i- 



(26) 



If we adopt our limiting case that Hm > 0.1 in order to 
be consistent with various measures of the contemporary 
matter density, we find that the restriction F < 0.46 from 
(H) asserts that 



MkevY< 1.2 X 10" 



(27) 



We illustrate this bound in Figure ^ by the vertical 
boundary for Tyr^ 10^. Lastly, because (1 -I- zeq) ~ 10^ 
in this scenario and fix is not more than a factor of five 
larger than our lower bound on flM , requiring jmin^ 10'^ 
provides only a weak restriction on X lifetimes. It may 
be more useful to take advantage of the large scale fea- 
ture that would be present in the power spectrum due to 
the injection of relativistic energy in order to limit relic 
properties. We do not explore such bounds as this would 
require specifying MkeV and Y separately. 
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C. SNIa Constraint on Relic Decays 

W ith t he SNIa bound on relativistic energy from Sec- 
tion |II C| , it is now easy to obtain a SNIa bound on the 
rehc decay properties. Using (15) and the 95% upper 
hmit in (o), we have 



M, 



keV 



Y^Tyr< 5.3 X 10^ 



(28) 



where we have once again assumed h < 0.88. Notice 
that this follows from a bound on the contemporary 
RED and, therefore, is most pertinent to late decays 
(i.e., Tyr^ 10'*). Again, this constraint is shown in Figure 
where our bounds on the properties of decaying relics 
are summarized. 



IV. CONCLUSIONS 

Constraints on the cosmological RED can provide a 
fundamental probe of particle physics beyond the stan- 
dard model. In this paper we have discussed constraints 
on the RED during four distinct epochs arising from 
BBN, the CMB, LSS and SNIa (see Figure 2). Fur- 
ther, we have shown how these bounds constrain the 
mass and lifetime of a hypothetical big bang relic (see 
Figure 5). Somewhat surprisingly, the RED at the cur- 
rent epoch is relatively unconstrained: we have shown 
that the magnitude-redshift relation for SNIa is consis- 
tent with a flat universe comprised of up to 20% rela- 
tivistic energy. Conventional wisdom suggests that the 
RED today must be small to allow for sufficient growth of 
large scale structure and not appreciably alter the CMB 
anisotropy power spectrum. The LSS bound does, in 
fact, significantly limit the RED near the epoch of matter- 
radiation equality. Any RED consistent with subsequent 



growth of LSS, redshifted to the current epoch, would be 
quite small. Conversely, a RED that is large yet accept- 
able with respect to the SNIa bound would clearly inhibit 
the growth of LSS if it were redshifted to the past. How- 
ever, a long-lived particle that decays sufhciently late as 
to avoid the LSS constraint could nevertheless contribute 
substantially to the RED today. In this case, the relevant 
constra int wo uld be the aforementioned CMB bound (see 
Section [II A). Avoiding the CMB constraint requires the 
X particles to be very long-lived if they are to contribute 
appreciably to the RED today. Consider a big bang relic 
with MkevY — 1-2 X 10^'^(a 30 eV neutrino is an ex- 
ample of a particle with the necessary abundance) and 
a very long lifetime, Tyr = 2 x 10^. The decay prod- 
ucts of this relic contribute fl^h^ ~ 0.1 and its proper- 
ties are marginally consistent with the growth of LSS . 
In addition, the decay products are produced sufficiently 
late so as to contribute an unconstrained ISW perturba- 
tion at low multipole moments, peaking around £ 10, 
and to change the angular diameter distance to the last- 
scattering surface by only ~ 7% (see Figures 3 and 4 of 
Kaplinghat et al., |^). These effects cannot be ruled 
out by current CMB data. In Figure || we show the evo- 
lution of the RED including the decay products of this 
hypothetical, long-lived big bang relic. 
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